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Abstract 



A set of world-line deviation equations is derived in the framework of Mathisson-Papapetrou-Dixon 
desciption of pseudo-classical spinning particles . They generalize the geodesic deviation equations. We 
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0" 1 1 Introduction 
u 

bJ[), The geodesic deviation equation , 

; -^-r = -R^wv^vP (i) 

H , 

is one of the well-known equations of general relativity. This equation shows the role of the space-time curvature 
on the motion of test particles and has important applications , namely, it is used for calculating relative accel- 
erations of nearby particles in an observer-independent manner, and may be integrated to give the Lyapunov 
exponent in the study of chaotic behaviour of particle's orbits. 

There are several ways to derive this equation. One standard derivation is as follows: Consider a two 
parameter family of geodesic x^(t\ A) in which r is the parameter along the orbits and A characterize different 
orbits. To obtain the separation of two nearby geodesies, one of them is taken as "fiducial" geodesic described 
by, say, x M , and the other by x^ + An M where An M represents the separation. The equation for ?i M is then 
obtained by inserting x^ + An M into the geodesic equation, comparing terms linear in A and neglecting 0(A 2 ). 
This idea have been used in to obtain generalized geodesic deviation equations by considering expansions 
containing higher orders of A. These generalized equations have been applied, for example, in [1 to the problem 
of closed orbital motion of test particles in the Kerr space-time and in [2] to the orbital motion in Schwarzchild 
metric. Other interesting generalizations may be found in [3j and 0]. Another way to derive the equation 
is by varying a suitable action [S]. This method was used in to generalize the equation to the so called 
" string deviation equation" . The quantization of the geodesic equation was discussed in 7\ and generalizing 
the equation to Kaluza-Klein theories was done in M. Some other aspects or applications of the equation have 
been discussed in 0. 

When forces other than gravity are present, or the particle has some internal structure, it would no longer 
move along geodesies in general. In these situations a "world- line deviation" equation may be obtained by 
modifying the geodesic deviation equation by taking the effect of the matter field or internal structure into 
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account. Thus using methods described in the previous paragraph, a world- line deviation equation was ob- 
tained in mil for the of motion of charged particles in the framework of Einstein-Maxwell theory. Another 
generalization of this type was made in (see also f° r describing spinning particles and was used there 
to study epicycles. 

The equations used in jll) to determine the particle's trajectories are a simplified version of Mathisson- 
Papapertou-Dixon (MPD) equations [T3]. The later equations have been widely used to study the motion 
of spinning particles (e.g. see |14j and references therein). The aim of the present work is to obtain a 
world-line deviation in the framework of MPD description of spinning particles. In the following sections 
we apply the prescription mentioned earlier to obtain the equations, and provide an example in which they 
could be integrated analytically. Throughout the work, 1, 2, 3, 4 would stand for u, v, x, y as indices , R^uap — 
dcX^vp — dpT^ ua + Y^ a $Y 5 v p — T^psT 6 ^ , and V represent covariant derivative, the space-time signature 
would be ( — h ++), and [fj,, u] stands for [w — vfj,. 

2 The equations 

The motion of a spinning particle is described by MPD equations 

n n c i 

-^T = -^ukxVs k \ (2) 

\- = (3) 
V^ u = (4) 

in which — is the particle's four-velocity, p^ its four-momentum and s Ml/ its spin tensor. We confine 
ourselves to time-like orbits with 



-1. (5) 



MPD equations guarantee that the particle's mass and spin are conserved 



p^p^ — constant = — m 2 , (6) 



s^vS^ = constant = s 2 . (7) 



Now consider a set x^r; A) describing the world- lines of spinning particles of the same spin-to- mass ratios and 
define 



DX 

and 



DX 

Inserting x M (r) + An Ai (r) into an d looking for A terms, we obtain the following equations 

ijjn i Dn v 1 1 

-=T- = -R»p aK v K n a pP - -R^ ap -— S ^ - -R» ua pv»J a P - -V K R" ua pn K v v (8) 
Lit 1 Lit Z i 

^— = s ^R^ KaP n^+p^^-+jK^ (9) 
Lit Lit 



W + -W = ( 10 ) 

respectively. If we turn off the spin, all terms except the first in the right hand side of JHJ vanish and the 
geodesic deviation equation results. If we set p^ — mv 11 equations (f%))- (|10|) reduce to those of [H]. It can be 
seen from the above equations that there is no evolution equation for it should be obtained indirectly from 
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equations (|8|)- (|10fl . The situation resembles the case of MPD equations in which no direct equation exists for 
The following equations 

« M -^=0, (11) 

= 0, (12) 

s^J^ = 0. (13) 

are helpful in this regard. They stem from 10-10 respectively. The latter two equations can be used for 
world-lines describing nearby particles of the same spins and masses. 

3 The motion in a GW space-time 

In this section we consider the world-line deviations of spinning particles in the space-time described by the 
following metric 

ds 2 = -dudv - K(u, x, y)du 2 + dx 2 + dy 2 (14) 

in which u = t — z,v — t + z are light-cone coordinates and K(u,x,y) = f(u)(x 2 — y 2 ) + 2g(u)xy. This 
metric represents a plane gravitational wave of arbitrary polarization and profile characterized by f(u),g(u) 
propagating in z-direction. In this space-time, the MPD equations admit the following solution 

v» = (1,1,0,0), 

p» = (m, m, 0,0), (15) 
s 1 " = s 2 ^ = 0, s 3i = <r. 

This solution describes a particle of mass m and of spin a (s M := ^= 6^ VK pP v ' s Kp , with e 1234 = —1) sitting 
in the origin of the coordinates with its spin along the z direction. For the above world-line, the geodesic 
deviation equation leads to the following equations 

=-f{u)n\r)~g(u)n\T) 1 (16) 

^P =/(«)n*(r)-5(«)n 8 (r). (17) 

We now consider two particles of the same spins and masses , one initially at the origin and the other at a 
nearby point with a specific separation. We aim to calculate this separation at any value of the parameter r. 
We take (|TK|) with (r, r, 0, 0) as the fiducial world-line. Now equation (fTTTf) results in 

J 12 (r) = 0, (18) 

J^(r) + J 23 (r) = -^/(r), (19) 
m 

J^{t) + J 2 \t) = ^f(r). (20) 



Equation (JT^-O results in 



dn 1 ^) dn 2 {r) 



= 0, (21) 



dr dr 

i X (r)+i 2 (r) = 0, (22) 
J 34 (r) = 0, (23) 
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respectively. From equation |jSJ) we obtain 

j (t) = constant = a, (24) 

j (r) = constant = (3, (25) 

: f(u) (J 13 (r) - mn 3 (r)) + g(u) (J 14 (r) - ran 4 (r)) , (26) 



dr 



dr 

Equation © leads to 

dJ 12 (r) (dn 2 (r) dn^r) 



f(u) ( J 14 (r) - mn\r)) + g(u) (j 13 (r) - mn 3 (r)) . (27) 



dr ™h^-^> +2a < 

^ = m^-fir), (29) 

dr dr 

dJ 14 (r) dn\r) , 4 

= W ^r-~ J(T) ' (30) 

= 2a 9 ( U )n 3 (r)-2 ( r/( U )n 4 (r)+m^^- J 3 (r), (31) 

= -2 ( T 5 ( M )n 4 (r)-2a/( W )n 3 (r)+m^^- J 4 (r), (32) 

dJ 34 (r) 



= 0. (33) 



Now, comparing l|24ll . 1251) . and (1221) we get /3 = —a. If we use the same parameter to describe both of the 
world-lines, then n* = 0. So we set n 1 (r) = — n 2 (r) which is consistent with (|21|l . Taking this into account, 
equations p8|) and 128|) give us — n (r) = n 2 (r) = — r + 7 in which 7 and a are the initial values of n 2 

and respectively. The transverse components of n M can be obtained from equations i|19|) - (|20[l . i|26|) . (|27|l . 
and l|29|) - (|33|) if one knows the second particle's world-line and spin orientation. Two interesting cases are as 
follow. 

Case 1: The second particle moves on a nearby geodesic. We set J 13 (r) = 0, J 14 (r) = 0. Thus j 3 (r) = 
J 4 (r) = and J 23 (r) = - 2( r^, J 24 (r) = 2a^. It follows that 

d 2 n 3 (r) 
d 2 n 4 (r) 

±±p. =f{u)n\T)-9{u)n\r), 

as we expected. 

Case 2: The second particle's spin is such that j 3 (r) = j 4 (r) = 0, and J 13 (r) = mn 3 (r), J 14 (r) = mn 4 (r). 
It follows that 

= ^(/(„) n V)-fl(«)» 3 (r)), (34) 
dr m 

^ = -(/(^W+jW^r)), (35) 
dr m 

which can be solved for n 3 (r),n (r) if f(u),g(u) is given explicitly in terms of w = r. 

4 Conclusions 

The equations we have found describe the world- line deviations in the framework of MPD equations determining 
the effect of the spin-curvature coupling on relative accelerations of nearby particles. They may be helpful in 
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different applications including the study of chaotic behaviour of spinning particles in certain space-times |15| . 
Another important application of these equations is to find approximate solutions to MPD equations by using 
a known solution, an example is the case two of the previous section. We applied the equations to the case of 
motion in the space-time of a gravitational wave and showed that they can be integrated analytically. A linear 
approximation of these equations may be useful in some applications. In more complex situations they may be 
solved at least numerically. In situations in which a higher accuracy is needed one can use more sophisticated 
equations containing higher orders of A in a systematic way. If the particles experience extra forces like the 
Lorentz force, the equations can be modified by adding suitable terms. The example of the previous section 
deserves a more extended study. 
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